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Abstract 



Ph! Neutrinoless double beta decay (/3/3)ov occurs through the magnetic coupling of di- 

mension five, A^/m„*, among the excited electron neutrino u*, electron and W boson if 
u* is a massive Majorana neutrino. If the coupling is not small, i.e., \ > 1 the mass of 
the excited neutrino must not be gless than the Z boson mass, mz- Since v* contributes 
in the (/30) Ol/ decay as a vertual state, this decay will give an oppotunity to explore the 
much heavier mass region of v* . 

In this paper, we present the decay formula of (PP)ou decay through the v* exchange 
and discuss the constraint on the coupling constant and the mass of the excited neutrino. 
By comparing the recent data for 76 Ge, we find A^ (lTeV/m^*))(mjv/lTeV)3 < 4.1-10 -3 
where is the Majorana mass of the excited electron neutrino. If mjv = rn v * and 
> 1, we find the mass bound for the excited Majorana neutrino as m v * > 5.9-10 4 TeV. 
In order to obtain the constraint on the composite scale A, we have to specify the 
model. For the mirror type and the homodoublet type models, A^* ] jm Vif = //(a/2A) 
*e-mail address: takasugi@phys.wani.osaka-u.ac.jp 



where / is the relative strength of gauge couplings. Then, we obtain A > 170 f(m N /neV) ^TeV. 
For the sequential type model, A/toj,* = fv/(\^2A 2 ) where v is the vacuum expec- 
tation value of the dopublet Higgs boson, i.e., v =250GeV. In this model, we find 
A > 6.6/^(m 7V /lTeV)3TeV. 
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1 Introduction 



If neutrinos are composite particles, there exist the excited neutrinos which couple to the 
ground state leptons by the dimension five magnetic coupling |],[?]. This interaction is 
expressed as|| 

\("*) 

L int = g^ea^( V * L R + ^L^d^W' + h.c, (1) 
tti v * 

where u* is a heavy excited electron neutrino, L — (1 — 7s)/2, R — (1 + 7s)/2, m Vif is 
the mass demension which is of order the mass of v*. This interaction is derived by the 
SU(2) x £7(l)gauge invariant form and parameters r\ L and rj R should satisfy the following 
constraints by the chirarity conservation 0, 

\Vl I 2 + I Vr \ 2 = 1, VlVr = 0. (2) 

The extensive search of u* have been made by many groups @] and found that m Vif > 
91GeV by assuming that Xz > 1 which is the coupling for v* — > eZ decay similarly 
defined to . So far, the severe mass bound in the region m v * > mz is not obtained. 

The purpose of this paper is to explore the mass renge m Vif > rriz by using the neu- 
trinoless double beta decay (/9/?)oi/ by assuming that v* is a massive Majorana neutrino. 
Then the (PP)ou decay occurs through u* exchange. Since u* enters as a virtual state, 
we can investigate heavy v* . Panella and SrivastavaH was the first to investigate (J3(3)q v 
decay, but unfortunately their formula is incorrect. Here we shall present the correct 
expression for the decay . 

We start from the interaction in Eq.(l) which leads to the effective four point inter- 
action between leptons and hadrons as 

L eff = -G eff ea^( V * L R + V R L)u*d^l + h.c, (3) 

where Jl is the hadronic current and 

G eff = 2G F 9 -^. (4) 

m v * 

I 
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2 Decay formula of neutrinoless double beta decay 

In the second order perturbation of the effective interaction in Eq.(3), the (PP)ou decay 
takes place and the S- matrix for this decay is given by 

iq(x-y) 



S = i Glff 



j dxdy J dq 



m N e 



q 2 - m 2 N 



-e(x)a^(7^(i 7 £ l i2 + i 7 i 2 L) e cr (i/)^r(4(x)jJ(i/)), (5) 



2(2tt) 4 

where e c is the charge conjugation of e, i.e., e c = Ce T . In the follwoing, we take the 
S-wave of electron wave function which is given by 

\ 



< | e(x) | p >= iJ S (e)e iex ; ip s (e) 



le + m 
2e 



Xs 



t+m As J 



F (Z,e), 



(6) 



where e is the energy of electron and F (Z, e) is the relativistic Coulomb factor defined in 
Eq.(3.1.25) in Ref.*. The S-wave function is independent of the space coordinate. With 
this wave function, we obtain 



: U eff 

/ 2(2tt) 4 



J dxdy J dq 



e -ig( X -y) mN 

q 2 - m 2 N 



< N f | T(Ji(x)4(y) | N t > x 



t^(e 1 ,e 2 ,q°,q)e^ x ° + ^ - {e 1 ^ e 2 ) 



(7) 



where 



t^(e 1 , e 2 , q°, q) = M^^a^R + ^^(e^q, - e 2 g, )(q P + e l9p0 ) (8) 

Now we perform the q° integration first and then do the x° and y° integration. We 
found 



Rh — 



G eff 



dxdy I dq—e^-^ 



X 



^2(2^)3 J - J ' cj 
^\<N f \ 4(x) \nxn\ 4(y) \ N t > 



t I/a (e 1 ,e 2 ,uj,q) - (e x <-»• e 2 ) 



, (9) 



uj + E n - Ei + e 1 

where a/21 in the denominator is the statistical factor. In this expression, q° = uj = 



'm 2 N + q 2 with q —\ q \ which comes from the pole at q° = uj of the neutrino propergator. 
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The contribution from the q° = —uj pole turns out to be equal to the one from q° = uj 
pole in the approximation to keep S-wave of electron wave function. 

So far our formulae is exact up to the S-wave approximation of electron wave function. 
Next, we make following approximations: (1) We take the closure approximation where 
E n is replaced by the average value < E n > so that the sum of the intermediate states 
can be taken. (2)ei and e 2 in the numerator is neglected, (q^ — 62^0) (Qp + e i9po) — QpQp- 
(3) The energy denominators uj+ < E n > -E* + e x = uj+ < E n > -(^±^) + is 
replaced by fiom e =< E n > Here we replaced Ei and Ef by their masses Mj 

and Mf. The approximation (2) and (3) are valid because e\ and e 2 are of order lMeV 
which is much smaller than uj = \Jm 2 N + q 2 , the average value of q which is of order 1/R 
with R beeing the nuclear radius and [Xom e . 

With these approximations, we get 

R fi = jSll [ dxdy [ dq^e^-* < ^ 1 I * > x 

1 >/2!2(27r) 3 J J uj uj + /i m e 

q^s{e 2 ){a^\ a^fR + VrL)1>% fa). (10) 

Now we use the identity {a^, a pa } = 2(g pp g ut7 - g^ g vp - ie pvpu ^) and the non- 
relativistic approximation of the hadronic current as 

= E Jl(n)6(S-f n ), Jt(n) = T+(g v g, + g^oiW), (11) 

n 

where r„ is the position of the n-th nucleon in the nucleus and F(q 2 ) is the form factor 
defined by 

^ 2 > = (TTfeW' (12) 

with the value itla = 0.85MeV. We obtain 

G 2 ff r e !?> " m 
Rfi = r- Y \ dq— r < N f I Jl{n)Jl{m) \ N t > x 

m N {m 2 N g^ - q^^sfa^fR + rfgL^fa), (13) 

where r nm = r n — r m . It should be noted that m 2 N g ua , q°q° = uj 2 and qiq 1 terms contribute 
only to + — > + transition, while q°q j = uq j term contributes to + — > 0" transition. 



The + — > 2 + transition does not occur because we keep the S-wave of electron wave 
function. 

In the following, we concentrate on the + — > + transition. Sine we consider a heavy 
composite neutrino is much greater than rriA, we expand uj and uj + iitfm e in the q 
integration in the power of l/m N . We find 

+ ) 2 -*n'ff m ]£- + <° k J^}F\q\ (14) 

Thus we obtain 



{(m„ - ^ m e + ^^)M GTiJ v + ^((^) 2 M^ - ^-M' GT - \m' t )}. (15) 
77^ m N g A 3 3 



where 



M gt ,n — 


<N f \ 


E 








M' F = 


<N f \ 


E 








M' GT = 


<N f 


! E 



(16) 



-(+) T (+) ( 

n m \ 



T (+) T (+) 



R -)F 4 (x A ) \N t >, (17) 



f'nm 

. R 



(18) 



=< Nf | ]T r( + VW{3(? n • r nm ){B m ■ fnm) - °n • <? m }(^)i ? 5 (^) | ^ >, (19) 

where xa = mAr n m and neutrino potentials are 

F N (x) = ^{3 + 3x + x 2 )e~ x , F 4 {x) = ^(3 + 3x - x 2 )e~ x , 

F 5 (x) = ge- (20) 

After taking the spin sum and performing the phase-space integration, we find the 
half-life of the transition of the neutrinoless double beta decay due to the heavy composite 
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neutrino as 

7^(0+ -0+) = 4 f^V f^V G 01 | - W± + (/i ° me)2 )M CT „ 

\ m„ y Vm e / m N niA 

+ T±\iLfM' F - -M' - -M'\ | 2 , (21) 

where Goi is the phase space factor defined in Eq. (3.5. 17a) in Ref.6. The above expression 
is obtained by taking | rji | 4 + | tjr | 4 = 1 which is valid for both chirality cases because 
we imposed the chirality conservation in Eq.(2). That is, we consider only two cases, 
(Vl,Vr) = 0-,0) or (0,1). 



3 Constraint from neutrinoless double beta decay 

In the following, we analyse the constraint on the coupling by using the experimental 
half-life limit of 76 Ge which is measured by the Heidelberg- Moscow collaboration ||, 

T(0+ -> 0+ : 76 Ge) > 5.6 ■ 10 2 V 90%c./.. (22) 

We derive the constraint on composite parameters from this data. By using the values of 
nuclear matrix elements obtained by Hirsch, Klapdor-Kleingrothaus and Kovalenko[7], 

M GTyN = 0.113, M' F = 3.06 ■ 10~ 3 , M' GT = -7.70 ■ 10~ 3 , M' T = -3.09 ■ 10~ 3 , (23) 

and the phase space factor Gqi = 6.4 • 10 _15 /yr given in Ref.6, we find the half-life of the 
+ — > + transition for 76 Ge is given by 

T-i = ( W\* | (2* _ ^£ + + 7 .2 ■ 10- 2 ^ | 2 9.1 • 10-^/yr. (24) 

\ m vif J m A m A m N m A m N 

By comparing this formula to the data, we find by taking, 

/W| 2 m*L _ /fome + W + 7 2 . 1Q _ 2 rru . L4 . 1Q _ 8 (25) 
\ m w / m,A rriA M N rriA m N 

Our formula in Eq.(21) is valif for the case where the excited neutrino is a Majorana 
neutrino which is much heavier than = 0.85GeV. In this case, the dominant term in 
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Eq.(25) is the one which is proportional to the excited neutrino mass because iitfm e 
is of order lOMeV. Thus the constraint becomes a simple form as 



4 Limit on the excited neutrino mass and the com- 
posite scale 

Firstly, we shall discuss about the excited neutrino mass. If we assume that = m„*, 
we obtain 



This is the most stringent bound so far. 

Secondly we discuss the bound on the composite scale. Since the relation between the 
coupling constant and the composite scale is model dependent, we have to confine some 
specific models ||. 

1. Sequential type model 

We consider the sequential type model where the excited leptons are assigned as 



where [v* R ] indicates that there are two cases, with and without u R . 

There are two kinds of gauge invariant interactions. If v R exists, the interaction is 
given by 




(26) 



m u * > 5.9 • 10 4 TeV (A£ > 1). 



(27) 




(28) 



1 



■h^igf^W^ + g'f'jB, u )4>u R + h.c, 



(29) 



2A 
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where II is the left-handed ground state lepton doublet and is the doublet Higgs 
boson. This interaction gives the one in Eq.(l) with 



("•) 

with { VL , m ) = {1,0), (30) 



where v is the vacuum expectation value of ^/2(p which is about 250GeV and v* R 
mediates the decay. 

If v* R does not exist, the other type interaction will be present 



Um = 7^e R ^a^(gf^W« u + g'f^B^Ll + h.c. (31) 



This interaction gives 

*w f v 



m. 



V2A 2 



with ( Vl , Vr ) = (0,1) (32) 



and v* L meadiates. 



For both cases, we find 



A>6.6/^(^^)3TeV. (33) 
v lTeV ; V ; 



2. Mirror type 

In this model, the multiplet of the excited leptons are 



r * l * r * 



( * 
1 V* 



R 



The gauge invariant interaction is 



which leads the coupling 



\(" ) f 



and u* R mediates. In this case, we find 



(34) 



L int = ^LCJ^igf^W^ + g'f^B, u )L R + h.c, (35) 



("•) 

with (r] L , m ) = (l,0) (36) 



A > 170/(^)^TeV. (37) 



3. Homodoublet type 

In this model, the multiplet of excited leptons are 



V 



(38) 



R 



In this case, the gauge invariant interaction is the same as Mirror case so that we 
have the same coupling in Eq.(*) and v* R mediates. Thus, we obtain 



A > 170/(^)^TeV. (39) 

In the above, we obtain rather strong constraints on the mass of u* and the composite 
scale. These constraints may be modified if there is a mixing for excited neutrinos. There 
are two types of mixing; the flavor mixing and the v* L — v* R mixing. If such mixings are 
present, the excited neutrino mass m N should be replaced by a effective mass as 

m N =>• < m„, >= J2( U ej or V£j)mj, (40) 

3 

as the one appeared in the neutrino mass contribution to the ordinary neutrinoless double 
beta decay. The mixing parameters are defined by 

"l = T, u v N iL, 4 = E^^, (41) 

3 3 

where N* is the mass-eigenstate Majorana neutrino with the mass rrij. The effective mass 
< m vif > in essence comes from the Majorana mass of v* L or v* R . Thus, if the main part 
of v* mass comes from the Dirac mass, then the constraint is not valid. In particular, if 
v* is a Dirac particle, there is no constraint from the neutrinoless double beta decay. 
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